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The Gutzwiller projection is shown to induce nontrivial correlation between the spinon and the 
holon in the slave Boson theory of the t — J model. We find this correlation is responsible some 
subtle differences between the slave Boson mean field theory and the Gutzwiller projected wave 
function(GWF), among which the particle-hole asymmetry in the quasiparticle weight calculated 
from the GWF is a particular example. 
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Superconductivity results from Bose condensation of 
charged particles. In the BCS theory of superconductiv- 
ity, Ferniionic electrons are paired into Bosonic Copper 
pairs which then condense into a superfluid. Soon af- 
ter the discovery of high temperature superconductors, 
Anderson proposed an exotic way toward superconduc- 
tivity in this class of materials. Anderson's proposal for 
superconductivity is to fractionalize the electron, rather 
than pair them up[l|, 12, b|. According to Anderson's 
argument, doping holes into the parent compounds of 
the high temperature superconductors, which are two di- 
mensional antiferromagnetic insulators, would generate 
a liquid-type spin state called RVB state. Anderson ar- 
gued that the excitations above the RVB state are frac- 
tionalized. More specifically, the spin and the charge 
quantum number of the electron are carried separately 
by two kinds of excitations, namely a spin-^ Fermionic 
excitation called spinon and a charge 1 Bosonic charged 
excitation called holon. In such a spin-charge separated 
system, the charged holon is liberated from the Fermionic 
statistics of the electron and is ready to condense into a 
superfluid. 

A problem with Anderson's proposal is that the pre- 
dicted Tc is too high|j,|5j. The holon supercurrent is not 
efficiently dissipated in this mechanism. In the BCS the- 
ory of superconductivity, the supercurrent is dissipated 
by thermally generated quasiparticle excitations. These 
quasiparticles, which are charged Fermionic excitations, 
form a normal fluid and cause dissipative response in ex- 
ternal electromagnetic(EM) fleld. On the other hand, 
in an ideal spin-charge separated system, the Fermionic 
spinon has been deprived of the charge to cause dissipa- 
tive response in an EM field, while the Bosonic excitation 
of the holon system is very inefficient in dissipating the 
supercurrent. 

One suggested solution for the above problem is to 
recombine the separated spinon and holon (ij, [3, l3l- By 
spinon - holon recombination, the Fermionic spinon exci- 
tation acquire the charge to cause dissipation in an EM 
field, or, the charged holon regain the Fermionic statis- 



tics to be transformed into a normal carrier outside the 
condensate. 

However, the origin and even the exact meaning of such 
recombination are still elusive. In this paper, we address 
the problem of spinon-holon recombination in the slave 
Boson theory of the t — J model. The slave Boson theory 
is the most direct way to realize the RVB idea. In the 
slave Boson theory, the constrained electron operator Ci^ 
is rewritten as fiabl, hi which fi^ denotes the Fermionic 
spinon and b\ denotes the Bosonic holon. To be a faithful 
representation of the original electron, the slave particles 
should be subjected to the local constraint Yl,a flafi,'^ + 
b% = 1. 

At the mean field level, the local constraint on the slave 
particle is relaxed to a global one. The spinon and the 
holon then move independently. Wen and Lee proposed 
that the Gaussian fluctuation around the mean field so- 
lution, which are gauge fluctuations, tends to recombine 
the spinon and the holonjg. However, the lack of a re- 
liable treatment of gauge fluctuation in 2+1 dimension 
make it hard to make any quantitative prediction. 

In this paper, we address the problem of spinon- 
holon recombination in the Gutzwiller projected wave 
functions (GWF). The GWF is derived from the slave 
Boson mean field state by filtering out its unphysical 
components with doubly occupied sites. Technically, 
the Gutzwiller projection, which enforces the local con- 
straint, accounts partially the effect of gauge fluctua- 
tion. Thus some kind of spinon-holon recombination is 
expected in GWF. We find the Gutzwiller projection not 
only enforces the local constraint, but also induces longer 
range correlation between the spinon and the holon. This 
correlation is given by a well defined correlation function 
and is found to be responsible for some subtle differences 
between the mean field theory and the GWF. For ex- 
ample, the particle-hole asymmetry in the quasiparticle 
weight calculated from GWF is shown to be due to this 
correlation. 



We start from the t — J model which is defined as 
H = -t ^ {c\„cj,a + h.c.) + JY^ (S.,Sj - -n^n. 



<i,j>,a 



<i,3> 



(1) 
in which Si = Y.ap^l.a'^alJCi^p and Ui = Y.a^\,a^i,a- 

The electron operator Ci^a in (1) is subjected to the con- 
straint of no boule occupancy 



E4<.^'.-^i- 



(2) 



In terms of the spinon and the holon operator, the t — J 
model readsQ 

<i,j>,(7 

+ "^ / , UiaJiPJjpJja ~ JiaJiaJjf^Jjp), 
<i,j>,a,l3 

in which fia and bi are now subjected to the constraint 

E4-/^- + ^l^* = l- (3) 

In the mean field theory, the interaction term is de- 
coupled by introducing RVB order parameter Xij ^ 

'Eaiflafj,'^) and Ay = J2ai<^c,fiflafj,p)^ i^ ^hich eaf} 
is the total antisymmetric tensor. The local constraint is 
relaxed to a global one. The mean field Hamiltonian for 
the spinon and the holon take the form 

H/ - E ^k/kV/k.^ + E ^k(/L/ik,^ + h.c.) (4) 



and 



H/,. = ^ ek6^&k 



(5) 



in momentum space. The mean field ground state is given 
by the following product 



G)mf = (6U)^1BCS}, 



(6) 



in which |BCS) denotes the BCS pairing state of the 
spinon. A mean field excitation, for, say, a hole-like 
quasipaticle, is given by 



|k,cr)MF = foJ^o7L|G)MF, 



(7) 



„t 



in which j^^ denotes the Bogliubov quasiparticle of the 
BCS Hamiltonian for spinon. 

The GWF for the ground state is given by project- 
ing |G)mf into the subspace that satisfy the constraint 
Eq.(3), 



\G)v,, = PG{bLo)^'\BCS), 



in which |G)var denotes the variational ground state 
and Pq denotes the Gutzwiller projection. Similarly, 
PG|k)MF provides a variational guess for the quasipat- 
icle excitation on |G)var- The same construction for the 
quasiparticle excitation is used in some recent works J9j- 

The spinon-holon recombination in GWF can be most 
easily seen by inspecting the quasiparticle weight Z^. In 
the mean field theory, Zk scales linearly with hole density 
, since the coherent spectral weight is caused by holon 
condensation. However, Zk calculated from GWF van- 
ishes more slowly near half filling. For example, Zk cal- 
culated from the one dimensional GWF vanishes as x^^"^ 
near half filling'lO'|. Similar behavior is also observed in 
two dimensional systems llj. In some cases, Zk calcu- 
lated from GWF can even be nonzero at half filling|l2l|. 
This discrepancy can be naturally explained by invoking 
spinon-holon recombination in GWF. 

Now we define spinon-holon recombination in GWF 
more explicitly. For clarity's sake, we first consider a 
half filled system and assume the spinon pairing term to 
be absent [13 ■ The mean field ground state of the system 
is then given by a half filled spinon Fermi sea, while a 
hole-like quasiparticle is generated as follow 



|k,T)MF = 6i=o/-k,i|FS), 



(9) 



in which |FS) denotes the half filled Fermi sea. At the 
mean field level, |k, t)MF is interpreted as a state with an 
added holon and a hole of spinon on the background of 
|FS). In real space. 






(10) 



Thus, the added holon and the hole of spinon move in- 
dependently on |FS} and have a chance of only -i- to 
recombine into the original electron. This explains the 
vanishing of the quasiparticle weight at half filling in the 
mean field theory. 

Now let's see what happens when the local constraint 
is enforced by the Gutzwiller projection. The projected 
wave function for the above quasiparticle is given by 



|k,T)va.. = ^E^*^" 



-^)PGfel/,i|FS). (11) 



(8) 



Following the convention of the mean field theory, we in- 
terpret the index i and j as the locations of the added 
holon and the hole of spinon. The difference between 
k, t)var and |k, Dmf can be understood as follows. In 
the mean field theory, the relative motion between the 
added holon and the hole of spinon is described by plane 
wave factor. The background particles in |FS) do not 
contribute to their correlation. After the Gutzwiller pro- 
jection, thees added particles get correlated with the 
background particles through the local constraint. Thus 
the background particles in |FS) can also contribute to 



the correlation between the added holon and the hole of 
spinon in |k, t)var- 

First, we assume i ^ j. Due to the local constraint, 
spinons in |FS) must keep away from the site i, which 
is already assigned to the added holon. At the same 
time, site j should be doubly occupied in |FS), since it 
accommodates a spin in the projected state. The sites 
other than i and j should be all singly occupied in |FS). 
Thus, the probability of finding the holon at site i while 
the hole of spinon at site j, which is denoted as P,;j, 
is given by that of finding site i empty, site j doubly 
occupied, and all other sites singly occupied in |FS). For 
J = J, in which case a real electron is removed from site 
i, the corresponding probability (Pij) is given by that of 
finding site i occupied by a down spin and all other sites 
singly occupied in |FS}. The ratio between Pij and Pa is 
given by 
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EJV'al 



(12) 



Here, Va denotes the amplitude of |FS) for a general 
configuration |a) in the subspace of no double occupancy, 
Tpfi denotes the amplitude for the configuration derived 
from I a) by moving an electron(of either spin) from site 
i to site j. This statistical sum can be evaluated easily 
with Variational Monte Carlo method. 

To see the value of the thus defined correlation func- 
tion, we consider some special cases in which the effect of 
spinon-holon recombination is manifest. As a trivial ex- 
ample, we consider the case of removing an electron from 
a fully polarized spin background. According to the slave 
Boson mean field theory, the quasipaticle weight is given 
by the hole density and vanishes at half filling. However, 
since the spin is fully polarized, the state in considera- 
tion is in fact a mean field state in which the quasiparticle 
weight should be exactly one. This discrepancy can be 
easily resolved by noting that Pij is nonzero only for i = j 
on a fully polarized spin background. This indicates that 
the added holon and the hole of spinon must occupy the 
same site and thus recombine into a real electron on that 
site. 

As a less trivial example, we consider a state with an- 
tiferromagnetic long range order. As we will show below, 
Pij decays exponentially in this case. The exponential de- 
cay of Pij can be understood by noting the opening of a 
SDW gap between configuration \(3) and |a) in Eq. (12). 
Thus the added holon and the hole of spinon will form 
a well defined bound state in such a spin background. 
Since the bound state has a nonzero overlap with a bare 
electron, the quasiparticle weight should be finite even at 
half filling. Calculation in 'l2| does find a finite quasi- 
particle weight for such a state. 

Now we present the numerical results. Figure 1 shows 
the correlation function for the one dimensional projected 
Fermi sea. The correlation function decays as l/r at large 
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FIG. 1: Spinon - holon correlation function for the one di- 
mensional projected Fermi sea. The inset show the data in 
logarithmic scale. 



distance and is unnormalizable. Thus the quasiparticle 
weight should vanish at half filling. The result for the 
two dimensional d-wave RVB state is shown in Figure 
2. The correlation function now decays as l/ra at large 
distance, which is also unnormalizable. Thus the quasi- 
particle weight should also vanish at half filling in this 
state. In Figure 3, we show the result for the antiferro- 
magetic ordered state. As mentioned above, the correla- 
tion function decays exponentially in this case. However, 
the short range behavior is quite similar to that of the 
d-wave RVB state. 
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FIG. 2: (a) Spinon - holon correlation function for a d-wave 
RVB state. The d-wave RVB order parameter A/i = 0.25. 
The calculation is done on a 20 x 20 lattice, (b) Power law 
decay of the spinon - holon correlation function. L denotes 
the linear size of the lattice, rmax is largest distance that can 
be defined on a periodic lattice of size L. 
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FIG. 3: (a)Sanie as Figure 2, but now for an antiferromagnetic 
ordered state. The SDW and d-wave order parameters are 
AAp/t = 0.1 and A/i = 0.25. (b)Exponential decay of the 
spinon - holon correlation function. 



The spinon- holon recombination discussed above refers 
to the correlation between the added particles rather 
than that between the background particles. At finite 
doping, a holon condensate is established and the back- 
ground spinon has a chance of x to behave as a coherent 
quasipaticle. In fact, the latter is the only way that the 
system can build a finite quasiparticle weight on the par- 
ticle side of the spectrum. Following essentially the same 
steps detailed above, one easily sees that the added slave 
paricles(or hole of them) during the process of injecting 
an electron have no correlation with each other. Since the 
quasiparticle weight due to holon condensation is by def- 
inition particle-hole symmetric, the particle-hole asym- 
metry in the quasiparticle weight should be attributed 
to the spinon- holon recombination effect discussed in this 
paper [ll|. 

Especially, while the quasiparticle weight for adding an 
electron is constrained by the sum rule to vanish at half 
filling, the quasiparticle weight for removing an electron 
can be finite at half filling, provided that the added holon 
and the hole of spion form well defined bound state, as is 
the case in the antiferromagnetic ordered state. In such a 
case the quasiparticle can contribute a substantial part of 
the tunneling asymmetry, although for a full understand- 
ing of the latter one also should take into account the 
contribution from incoherent spectral weight, 14lj. When 
the added holon and the hole of spinon are less tightly 
bounded, the particle-hole asymmetry in the quasiparti- 
cle weight should be less dramatic, as is found in a recent 
work on the tunneling asymmetry of high temperature 
superconductors [m . 

We now discuss some other effects of spinon-holon re- 
combination. By recombining with a spinon excitation, a 



hole changes the local spin environment around it. This 
effect can be invoked to release the kinetic energy of the 
holes in certain spin background through spontaneous 
generation of spinon excitations in the system. An ex- 
ample in this respect is provided by a recent variational 
study on the shghtly doped cupratesp^. In that work, a 
t — J model with next-nearest-neighboring hoping term 
is considered. At small doping , it is found that a varia- 
tional state with spontaneously generated spinon excita- 
tions give lower energy than a state without spinon exci- 
tations. It is also found that the presence of these spinon 
excitations changes the charge correlation dramatically. 
This can also be understood in terms of spinon-holon re- 
combination defined in this paper. A detailed analysis of 
these issues will appear in separate paper. 

Finally, we revisit the problem of supercurrent dissipa- 
tion. In the presence of spinon excitations, a holon tend 
to bind with them and form a Fermionic composite ob- 
ject. A holon is thus transformed into a normal charge 
carrier out of its condensate. The superconductivity is 
gone when all holons in the condensate are transformed 
into normal carriers. At low doping, a small number of 
spinon excitation is enough to kill the superconductiv- 
ity which is not expected to alter the spin correlation 
significantly. This is proposed as a mechanism for the 
pseudogap in underdoped cupratesli||7|. 

In conclusion, an explicit definition for spinon-holon 
recombination is given on the GWF. The thus defined 
spinon-holon recombination is shown to be responsi- 
ble for the particle-hole asymmetry in the quasiparticle 
weight calculated from the GWF and some other subtle 
differences between the mean field theory and the GWF. 
A new mechanism for the dissipation of supercurrent in 
the RVB state is proposed based on the spinon-holon re- 
combination defined in this paper. 
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